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We  study  an  approximate  solution  of  the  Boltzmann  Equation  for  the  distribution  of 
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in  the  electric  field,  and  mean  energy  loss  raTe“To  collisions  are  discussed.  The 
importance  of  spread  of  energy  gain  and  loss  about  the  mean  rates  is  emphasized, 
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mean  momentum  transfer  collision  frequency,  but  on  any  time  scale  relative  to  the 
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ipThe  steady  state  spectrum  can  be  obtained  in  closed  form.  It  reduces  to  quadratures 
in  terms  of  the  momentum  transfer  cross  section  and  two  sums  over  energy  transfer 
cross  sections.  *  Its  connection  with  ordinary  diffusion  theory  is  pointed  out.  It 
explains  the  shortcomings  of  the  Continuous  Slowing  Down  Approximation. 


It..is  demonstrated  expLici tly  -for  rotational  excitations  that,  when  E  =  0,  the 
diffusion -approximatfon  reduces  to  the  correct  Maxwellian  distribution  at  the  gas 
temperature  in  the  limit  that  rotational  energy  spacings  become  small.  The  theory 
is  applied  to  Nitrogen,  where  energy  spectra  and  transport  coefficients  are  computed, 
and  compared  with  CSDA  calculations,  data,  and  more  accurate  calculations. 
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SECTION  1 
INTRODUCTION 


The  behavior  of  the  electron  energy  spectrum  in  a  weakly  ionized 
gas  in  the  presence  of  electric  fields,  momentum  transfer  collisions,  and 
energy  transfer  collisions  with  the  host  gas  is  a  long  standing  problem  in 
plasma  physics.  Historical  surveys  appear  in  References  1  and  2.  We  have 
in  mind  electrons  with  energy  from  zero  to  several  eV  or  a  few  tens  of  eV. 

The  electron  energy  spectrum  controls  many  processes  of  inter¬ 
est,  such  as  rates  of  chemical  reactions  (attachment,  detachment,  recom¬ 
bination,  etc.),  avalanche  rate,  electrical  conductivity  and  mobility, 
diffusion,  mean  electron  energy,  mean  momentum-  and  energy-transfer  col¬ 
lision  frequencies,  etc.  Weakly  ionized  gases  occur  frequently  in  natural 
and  laboratory  situations:  for  example  the  ambient  or  nuclear-disturbed 
ionosphere,  and  in  ionization  following  dosed  air  due  to  gamma-rays  (e.g. 
nuclear  EMP),  or  X-rays,  or  the  passage  of  fast  charged  particles.  If  an 
electric  field  is  present  it  can  cause  the  electron  cloud  to  drift,  heat 
it,  and  possibly  produce  avalanching. 

HISTORICAL 

The  spectrum  is  controlled,  of  course,  by  the  Boltzmann  Equa¬ 
tion,  and  techniques  for  solving  it,  and  certain  progress,  were  establish¬ 
ed  decades  ago  (e.g.,  References  3  and  4).  Useful  answers  were  hampered 
by  lack  of  basic  molecular  cross  sections,  and  the  difficulty  of  analytic¬ 
ally  solving  the  Boltzmann  Equation.  An  expansion  in  Legendre  Polynomials 
in  the  angle  between  the  velocity  vector  and  the  electric  field  is  useful 
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in  the  swarm  regime,  in  which  the  drift  velocity  is  small  compared  with 
thermal  velocity,  and  this  expansion  will  be  employed  in  the  present  work, 
at  least  in  its  early  stages. 

Harder  to  handle  is  the  collision  integral,  because  of  discrete 
energy  loss  in  inelastic  collisions.  The  integral  can  be  expanded  in 
powers  of  mean  energy  loss  per  collision  divided  by  the  average  swarm 
energy.  When  the  expansion  is  terminated  with  the  first  non-trivial  term, 
representing  mean  energy  loss,  the  ensuinq  results  are  called  the  Contin¬ 
uous  Slowing  Down  Approximation  (CSDA).  It  involves  an  energy-weighted 
sum  over  excitation  cross-sections  called  the  "Energy  Loss  Function",  or 
the  "Stopping  Cross  Section",  measuring  the  mean  enerqy  loss  rate  of  mono- 
energetic  electrons. 

Data  available  more  recently  have  made  it  possible  to  construct 
the  energy  loss  function  at  these  relatively  low  energies  where  it  is  ex¬ 
tremely  difficult  to  calculate  from  first  principles.  (At  higher  ener¬ 
gies,  >  1  keV,  the  Loss  Function  can  be  calculated  in  the  Born  approxima¬ 
tion,  where  it  is  known  as  the  Bethe  mean  stopping  power  formula.)  Using 
it,  energy  spectra  and  bulk  plasma  parameters  computed  in  the  CSDA  have 
been  shown  to  disagree  with  experimental  values  by,  for  example,  perhaps 
30%  for  mobilities,  and  by  large  factors  for  the  avalanche  rate  at  low 
temperatures  (e.g..  Reference  5).  Poor  agreement  with  avalanche  indicates 
that  the  CSDA  is  particularly  deficient  in  reproducing  the  high  energy 
tail. 

In  the  special  case  of  constant  mean  free  path  and  energy  loss 
only  to  elastic  recoil  (applicable  for  example  in  noble  gases  at  energies 
below  electronic  excitation),  the  transport  coefficients  computed  in  the 
CSDA  agree  reasonably  well  with  experiments.  The  spectrum  in  this  case  is 
known  as  the  Druyvesteyn  spectrum.  Historically  it  was  the  first  non- 
Maxwellian  spectrum  computed  for  a  plasma  in  the  presence  of  an  electric 


With  the  advent  of  better  experimental  cross-sections  and  com¬ 
puting  facilities  attention  turned  to  numerical  solutions  of  the  complete 
Boltzmann  Equation  with  inelastic  collision  integral,  by  either  Monte 
Carlo  or  finite  difference  techniques.  This  enabled  relatively  accurate 
energy  spectra  and  swarm  parameters  to  be  computed,  which  compare  general¬ 
ly  quite  successfully  with  experiments.  Indeed,  one  way  of  determining 
basic  electron-molecule  cross-sections  is  to  assume  a  form  for  them 
(guided  by  previous  experiments  and  quantum-mechanical  calculations),  and 
compute  the  resulting  swarm  parameters.  The  cross-sections  are  then 
adjusted  until  computed  parameters  agree  with  experimental  ones.  This 
technique  has  been  successfully  employed  by,  for  example,  A.  V.  Phelps  and 
co-workers  over  some  years  (for  example.  References  6  and  7). 

Thus  numerical  methods  exist  to  compute  spectra  and  transport 
coefficients.  These  numerical  techniques  do  not  offer  physical  insight 
into  swarm  behavior,  and  often  do  not  suggest  physically-motivated 
approximation  schemes  that  might  be  useful  in  specific  applications. 

APPROXIMATE  APPROACH 


Here  we  prefer  to  study  more  intuitively  the  processes  that 
determine  the  spectrum.  Our  goal  is  not  precise  solutions  of  the  Boltz¬ 
mann  Equation.  Rather,  we  build  on  physically  motivated  models  to  con¬ 
struct  an  approximation  to  the  inelastic  collision  integral  superior  to 
the  CSDA.  This  will  lead  us  to  an  approximate  differential  equation  for 
the  time  evolution  of  the  spectrum  in  a  time  varying  field,  and  to  a 
closed  form  analytic  expression  for  the  steady  state  energy  spectrum  in 
terms  of  the  host  gas  momentum-transfer  cross-section  and  two  weighted 
sums  over  inelastic  energy-transfer  cross-sections.  The  availability  of 
basic  electron-molecule  cross-sections  in  interesting  gases  makes  the 
approach  advantageous  and  useful. 


The  argument  makes  clear  the  importance  of  spread  of  energy  loss 
about  the  mean  (energy  diffusion,  loosely  speaking).  Usual  treatments, 
such  as  the  CSDA,  which  retain  only  the  mean  energy  loss  rate,  miss  a  cri¬ 
tical  piece  of  the  physics  determining  the  spectrum.  The  approximate 
soectrum  retaining  this  critical  physics  can  be  reduced  to  quadratures, 
and  is  essentially  no  more  complicated  in  form  than  the  CSDA.  It  makes 
clear  the  physical  reason  for  the  inadequacy  of  the  CSDA,  and,  in  particu¬ 
lar,  it  becomes  clear  why  the  tail  is  so  low  in  that  approximation.  We 
find  the  transport  coefficients  in  Nitrogen  come  out  in  closer  agreement 
with  experiment  than  in  the  CSDA.  Other  gases  will  be  studied  in  the 
future.  Hopefully  the  accuracy  of  the  computed  coefficients  in  these 
gases  will  be  sufficient  for  a  number  of  applications. 

The  validity  of  our  approximation  (which,  as  a  differential 
equation,  is  local  in  energy  space)  is  still  limited  by  discrete  finite 
energy  loss  jumps  in  inelastic  processes  such  as  electronic  excitations. 
We  require  energy  steps  to  be  small  compared  to  the  spectrum  width. 
Examples  of  its  limitations  will  be  seen  in  the  avalanche  rate  where  our 
imperfect  treatment  of  energy  loss,  although  improved  over  the  CSDA's, 
results  in  a  somewhat  high  rate,  and  in  the  N2  drift  velocity  at  a  few  eV 
where  the  tail  of  the  spectrum  can  excite  electronic  states,  suddenly 
decreasing  the  electron  energy  by  6  or  7  eV  or  more.  However  vibrational 
excitations  cause  fractionally  small  energy  transfer,  and  the  approxima¬ 
tion  appears  to  be  excellent  when  vibrational  excitations  dominate,  or  at 
higher  energies  where  even  electronic  levels  are  not  much  more  energetic 
than  the  energy  spread  in  the  spectrum  itself. 

Aside  from  these  problems,  the  spectrum  we  deduce  reproduces  ex¬ 
perimental  drift  velocities  and  diffusion  coefficients  very  well  over  some 
5  orders  of  magnitude  in  the  experimental  parameter  E/p  (the  reduced  elec¬ 
tric  field)  spannina,  for  example,  essentially  the  entire  range  of  data 
compiled  by  Dutton  (Ref.  8). 


SECTION  2 

THE  BOLTZMANN  EQUATION 


The  Boltzmann  Equation  for  the  electron  distribution  function 
f(v)  in  an  electric  field  is 

-> 

f  =  C  +  Sources  -  Sinks  (1) 

at  m  v 

dropping  the  spatial  gradient  term  and  any  magnetic  field  from  the  left 
hand  side.  Here  e  is  the  magnitude  of  the  electron  charge.  Sources  may 
be  X-ray  ionization,  or  ionization  by  a  fast  charged  particle.  Sinks 
would  include  recombination  or  attachment.  Actual  sources  and  sinks  are 
problem  specific,  and  are  not  considered  further  here.  We  concentrate  on 
the  effects  of  the  field  term  and  the  collision  integral  C  in  determining 
f . 


Many  bulk  plasma  parameters  depend  on  only  the  speed  distribu¬ 
tion  function 


F(v)  = 


(2) 


or,  equivalently,  on  the  energy  distribution  function  g(w), 

g(w)dw  =  F(v)v2dv, 
g(w)  =  ( v/m) F( v) 


(3) 


where  w  =  mv2/2  is  the  electron  energy.  If  n  is  the  number  density  of 
particles,  the  normalizations  are 


The  dimensions  of  g  are,  say,  cm-  eV~  \  The  equilibrium  Maxwell  distri¬ 
butions  for  a  gas  of  temperature  T  (in  energy  units),  which  obtain  when 
E  =  0,  are 


f M ( v )  =  n  r_[EL_'j 3/2  e-m2/2T 
M  2  t:T 


M(v)  =  4nf 

(w)  =  n  'w',1//2  p"w/T 

y  '  TtT  T  T1  6 


(5) 


LEGENDRE  EXPANSION  OF  THE  DISTRIBUTION  FUNCTION 


For  reduced  electric  fields  not  too  large,  the  departure  from 
isotropy  will  be  small,  and  f  is  usefully  developed  in  the  usual  Legendre 

-V  -V 

Polynomial  series  in  the  angle  9  between  v  and  E  (ref.  4) 


f  =  7  f  ,(v)  P  ,(cos9) 

=  f o  +  f l  cose  +  small  terms 


(6) 


The  quantitative  criterion 
collisions  dominate,  C  is 
series  are  inserted  in  Eq. 


is  V  .  <<  V.u  ,  .  When  el 

drift  thermal 

linear  in  f.  C  is  also  expanded, 
(1)  the  P ;  and  Pi  parts  become 


ectron-neutral 
and  when  the 


f0  controls  the  eneray  spectrum,  for  F  and  q  are 


F  (  V )  =  f  (fn+f^OSeJd^  =  4-rrf  0 

(8) 

a  =  ^lv  f0 
m 

while  fj  controls  the  bulk  electron  motion  with  a  drift  velocity 

v  =  I  r  f  v  d3v  =  1  f  (f 0+f ,cose)vd3v 
a  n  n 

(9) 

=  4tl  e  r  f  v3dv 
1 n  * 


LEGENDRE  EXPANSION  OF  COLLISION  INTEGRAL 

Thus  C0  is  the  contribution  of  the  collision  integral  to  the 
time  rate  of  change  of  g,  while  Cl  is  the  contribution  of  C  to  the  rate  of 
chanqe  of  that  part  of  the  spectrum  that  controls  the  swarm  momentum.  In 
qeneral,  for  a  weakly  ionized  qas  in  which  electron-neutral  collisions 
dominate  (ref.  4) 

ci  =  *  Vi  UO) 

where  N0v~m  's  the  '"ono-enerqetic  momentum  transfer  collision  frequen¬ 
cy,  is  the  momentum  transfer  cross-section,  and  N0  is  the  background 
aas  number  density. 

C-  receives  contributions  from  elastic  and  inelastic  processes. 

r  -  rplas  .  rinelas 


M 

T 


The  elastic  contributions  are  from  electron  energy  loss  to  heavy 
molecule  recoil  in  an  elastic  collision,  and  electron  heating  by  elastic 
collisions  with  the  host  gas  at  temperature  T  .  These  parts  are  (ref.  4) 

ce1as  =  i  _j__  2m  v  v^fQ  -t-  ?Ifl  (12) 

v  a(v2)  m0  m  m  a(v2) 

where  M0  is  the  molecule  mass,  assumed  much  qreater  than  the  electron  mass 
m.  The  first  term  in  Eq.  (12)  is  t*"  mean  energy  loss  rate  by  elastic 
collisions.  The  second  term  represents  the  heating  rate  by  the  warm  gas 
(the  Davydov  term).  Note  the  first  is  essentially  proportional  to  af0/aw, 
while  the  second  is  proportional  to  a2f0/aw2.  The  first  term  is  the  CSDA 
to  elastic  energy  loss,  and  the  second  is  a  diffusion-like  correction  due 
to  gas  motion.  Hiqher  order  corrections  (in  m/M0)  are  safely  negligible, 
and  the  CSDA  is  adequate  for  elastic  collisions. 

cjnelas  will  be  treated  in  the  next  section. 

SLOWLY  VARYING  FIELDS 

Eq.  (10)  is  inserted  into  Eq.  (7b), 


3f  2  eE  3  f  o  _ 

at  m  av 


(13) 


When  a/at  «  v  ,  we  have 

f l  =  eE_  3fo 

rnvm  a  v 
m 


(14) 


expressinq  the  first  order  distribution  function  in  terms  of  the  energy 
spectrum.  vm  is  greater  than  some  10  a  sec-1  in  sea  level  air,  so  condi¬ 
tion  (14)  holds  in  many  circumstances. 


To  obtain  the  equation  for  fo  we  combine  the  two  terms  in  the 
first  of  Eqs.  (7)  (2f i/v+3f i/dv  =  v" 2  d(v2fi)/}v),  and  insert  Eq.  (14)  for 
f  l : 


.  eV  X  fvi  £o,  .  c!'as  ♦  cr'as 

5t  3m  ;:v 2  ^  ^ 


(15) 


We  may  now  transform  to  an  equation  for  g(w): 

f  " 


23  +  1_  f e 2E 2  Cg  .  2  _d_  (Wg)  n  =  celas  +  Cine1as  (16) 

^t  "iw  m  v  3  ^w  g  g 

m 


In  this  equation  we  have  split  the  electric  field  term  into  two  parts,  one 
essentially  proportional  to  ^g/ ^w,  and  one  to  ^2gAw2.  We  also  define  the 
collision  integral  terms  on  g  rather  than  fo: 


^elas  _  2m  5 
g  Tw 


w[g  +  *  Tg  4 


7w 


)li 


(17) 


rinel as 
9 


is  taken  up  in  the  next  section. 


ENERGY  DIFFUSION  IN  THE  ELECTRIC  FIELD 


The  two  parts  of  the  electric  field  term  in  Eq.  (16)  correspond 
mathematically  to  a  mean  flux  (~dg/dw)  and  a  diffusive  flux  (~^2g/ 
3w2) .  That  is,  elastic  scattering  in  the  presence  of  an  electric  field 

produces  both  a  mean  gain  in  energy  (rate  ~  e  E  /mv  )  and  a  diffusion 


nvmvTw  u«uv<wv«iiv  w* 


9  9 

in  energy  (diffusion  coefficient  "  (2/3)we  E  /mvm).  That  elastic  (or 
momentum-transfer)  scattering  causes  velocity  diffusion  is  well-known. 

That  energy  diffusion  is  important  can  be  seen  by  comparing  its 
effect  in  Eq.  (16)  to  that  of  mean  energy  gain  on  a  Maxwellian 
distribution  of  temperature  T: 

-  1.  A_  ( wg ) 

Diffusive  Gain  =  3  5w _  =  +  2_  w  (18) 

Mean  Gain  g  3  T 

This  is  of  order  unity,  and  changes  sign  over  the  average  energy  3T/2. 
Diffusion  allows  particles  of  energy  w  >  3T/2  to  gain  energy  faster  than 
the  mean,  and  particles  with  w  <  3T/2  to  gain  energy  slower  than  mean, 
thus  spreading  out  the  spectrum.  Diffusion  adds  body  to  the  spectrum. 
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SECTION  3 

THE  INELASTIC  COLLISION  INTEGRAL 


The  general  form  for  the  collision  integral  is 
C  =  /  g ( w 1 )P( w' ,w)dw'  -  g(w)  /  P(w,w')dw' 


(19) 

l 

I 


where  P(w,w')dw'  is  the  probability  per  second  that  an  electron  with 
energy  w  collides  and  ends  up  with  energy  (w'.w'+dw1).  The  CSDA  to  the 
elastic  part  of  was  given  in  Eq(  1 7 ) .  For  the  inelastic  part,  colli- 
sional  excitations  to  discrete  molecular  states  k  of  energy  w^  and  cross 
section  c^(w),  contribute  to  P  an  amount 


P( w,w' )dw' 


=  No  v  l  <v 
k  * 


6(w-w^-w' )dw' 


(20) 


assuming  all  host  molecules  in  the  ground  state.  For  a  diatomic  molecule 
there  will  be  rotational,  vibrational,  and  electronic  states.  Excited 
rotational  levels  will  be  populated  thermally.  This  is  accounted  for  in 
Appendix  A  in  connection  with  the  evaluation  of  the  functions  L-n  and  M..^ 
of  Equations  (23).  For  simplicity  in  the  present  discussion,  we  write  N0 
as  if  all  molecules  were  in  their  ground  state.  The  integrals  in  Eq.  (19) 
are  then  readily  evaluated  for  these  inelastic  processes: 


This  clearly  exhibits  the  discrete  nature  of  energy  loss,  g  being  popu¬ 
lated  at  w  by  electrons  of  energy  w  +  w, .  Since  the  cross  sections  a 

k  • _ t  k 


are  not  exact  step  functions,  the  sum  C 


inel as 


itself  is,  of  course,  a 


continuous  function  of  w.  Developing  the  summand  in  a  Taylor  Series  about 
w,  we  exhibit  the  first  two  terms. 


.inel  as 


=  No  ~  rvgL,  +  -2_  (vgM.  )  1  + 

nW  1  1 '  i  1  >  < 


where 


L  .  =  y  w.  a,  (w) 

in  /-  k  k 


(23a) 


Min  =  U" k  °k(w) 


(23b) 


Lir}  is  the  Energy  Loss  Function  describing  the  mean  energy  loss  rate  to 
inelastic  processes.  According  to  Eq  (17),  the  elastic  part  of  the  loss 
function,  for  the  target  molecule  at  rest,  is 


L  i  =  iE  a  w 
el  m  m 


where  we  have  used  v  =  Nnva  .  The  total  mean  loss  rate  to  an  electron 

m  m 


of  energy  w  is 


L  =  L  ,  +  L. 
el  in 


Since  energy  loss  collisions  are  probablistic  processes,  there 
will  actually  be  a  distribution  in  the  energy  loss  about  the  mean  rate  L. 
This  is  approximated  by  the  term,  which  can  be  called  the  Stragqlinq 
Function  in  analogy  to  the  high  energy  case.  This  second  term  in  Eq. 
(22),  proportional  to  52/^w2,  approximates  the  true  discontinuous  random 


walk  in  energy  space  by  a  diffusion-like  process.  Contributions  to  the 
Straggling  Function  from  elastic  scattering,  M^,  are  smaller  than  by 

~  (m/Mo)w,  and  can  be  neglected.  M-n  therefore  essentially  does  not  dif¬ 
fer  from  the  total  Straggling  Function  M  =  +  M. 

The  CSDA  neglects  the  M..  term  and  retains  only  the  mean  loss 
rate  Lin  in  Eq.  (22) . 

Momentarily  neglecting  T  ,  Eq.  (16)  for  dg/at  now  reads 

l£  +  1_  {ifl!  [g  -  11-  (wg)  ]}  =  N0  ±_  [vgL  +  1_  (vgM)  ]  (26) 

bt  5w  m  v„  3  Sw  5w  "  5w 


or,  in  flux  conservative  form, 

13  +  12  =0 

at  aw 


where 


s  =  ®JL  f  g  -  |  1_  ( wq ) 

m  v  3  aw 

m 


-  NovgL  -  No  —  (vgM) 

aw 


is  the  flux  in  energy  space. 

MAGNITUDE  OF  SPREAD  IN  ENERGY  LOSS 

Like  L,  M  is  constructed  from  individual  excitation  cross 
sections.  The  ratio  M/L  is  likely  much  more  slowly  varying  than  L  or  M 
itself,  since  irregular  behavior  in  the  cross  sections  such  as  the  N2 
resonance  spike  near  2.5  eV  will  occur  in  both  L  and  M  and  divide  out  in 
the  ratio. 


Spread  about  the  mean  energy  loss  L  can  be  compared  with  L. 
From  Eq.  (28)  it  is  of  order 


Diffusive  Loss  -  1  3  (vgM)  -  M_ 

Mean  Loss  vgL  3w  "wL 


where  ¥  is  the  energy  over  which  Mg  =  (M/L)Lg  changes  much,  being  of 
order  the  swarm  temperature  or  a  character istic  energy  over  which  L  itself 
changes.  This  ratio  is  not  negligibly  small,  and  can  be  of  order  unity, 
showing  that  diffusive  contributions  to  energy  loss  can  be  comparable  to 
the  mean  loss,  as  was  the  case  for  energy  gain.  Since-  in  steady  state 
mean  loss  essentially  balances  mean  gain,  diffusive  loss  can  be  as 

important  as  diffusive  gain  in  determining  the  spectrum. 

IMPORTANCE  OF  ENERGY  DIFFUSION 

The  steady  state  spectrum  from  Eqs.  (27)  and  (28)  is  that  for 

which  $  =  0,  and  results  from  a  balance  between  total  energy  gain  from  the 

electric  field  and  total  energy  loss  to  collisions. 

It  is  instructive  to  temporarily  neglect  the  diffusive  contribu¬ 
tions  to  the  flux;  that  is,  the  M  term  in  the  collision  integral,  and  the 
term  -(2/3 ) 3 (wg) /3w  in  the  electric  field  expression.  The  steady  state 
spectrun  one  then  obtains  is  that  resulting  from  mean  energy  gain  e2E2/ 
mvm  balanced  by  mean  energy  loss  N0vL.  See  Figure  la.  This  spectrum 

reduces  to  a  delta  function  at  a  single  energy  Wj  where  the  rates  are 
equal.  Since  the  mean  gain  is  a  decreasing  function  of  w,  while  the  mean 
loss  is  an  increasing  function,  all  electrons  with  w  >  Wj  lose  energy 
faster  than  they  gain  it,  and  all  electrons  with  w  <  Wj  gain  energy  faster 
than  they  lose  it.  All  electrons  thereby  end  up  with  energy  wj.  Thus  the 
physical  process  which  actually  determines  the  energy  spectrin  is  spread 
in  energy  gain  or  loss  about  the  mean.  Mean  gain  and  loss  rates  determine 
only  the  average  energy  in  the  swarm. 
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MEAN  GAIN  RATE 


MEAN  LOSS  RATE 


RESULTANT  SPECTRUM 

Figure  la 

Bal ance  of  mean  gain 

and  loss  rates. 

Resultant  spectrun  is 

a  6  function  at  w 

COMPLETE  N 

/  ^COMPLETE 

/ 

REALISTIC 

SPECTRUM 


Figure  lb  Balance  of  complete  gain  and  loss  rates, 
CSDA  neglects  spread  in  energy  loss. 


^ 


It  is  clear  that  for  electrons  with  a  spread  in  energies  to  be 
in  steady  state  there  must  be  a  balance  between  gain  and  loss  rates  over  a 
spread  in  energies,  and  this  spread  determines  the  spectrum. 


This  shows  that  in  the  CSDA,  the  only  physical  process  which 
establishes  the  spectrum,  that  is,  gives  body  to  the  distribution 
function,  is  energy  diffusion  due  to  elastic  scattering  in  the  electric 
field,  since  in  the  collision  integral  only  the  mean  loss  rate  is 

retained.  Figure  lb  illustrates  the  CSDA  and  a  realistic  spectrum. 

Since  energy  diffusion  is  so  important,  we  must  keep  the 

diffusive  contribution  in  as  well  to  deduce  a  realistic  spectrum. 
While  L  is  the  lowest  order  term  in  the  expansion  of  C^,  it  affects  only 
the  average  swarm  energy.  M  is  the  lowest  order  contributor  from  the  col¬ 
lision  integral  to  the  energy  distribution. 

As  Figure  lb  shows,  the  CSDA  is  poor  because  it  retains  energy 
diffusion  in  the  E-field  energy  gain  term,  but  neglects  it  in  energy  loss 
processes.  That  is,  it  neglects  "half"  of  the  physics  determining  the 

spectrum. 

From  this  point  of  view  we  see  that  M  should  not  be  considered 
just  another  term  in  a  Taylor  Series  expansion  of  a  function.  Rather  it 
represents  a  critical  piece  of  physics  that  determines  the  distribution 
function,  and  must  not  be  neglected.  The  functional  form  in  which  M 
appears  is  a  Fokker-Pl anck  approximation  to  the  true  random  walk.  The 
approximation  will  be  best  when  the  energy  loss  per  collision  is  small 
compared  with  the  energy  spread  in  the  entire  spectrum.  Similar  to  the 

usual  Fokker-Pl anck  equation,  it  assumes  large  energy  changes  occur  as  a 
succession  of  small  ones. 


The  L,  M  approximation  to  Cg  is,  of  course,  local  in  energy 
space,  as  is  any  differential  approximation  to  a  random  walk  problem. 


n 

:i 


If  the  energy  loss  is  extremely  small  compared  with  w,  then 
Equations  (23)  show  M  «  wL,  and  the  CSOA  will  be  adequate.  This  happens, 
for  example  in  elastic  collision  recoil  as  mentioned  after  Eq.  (25).  If 
wk  is  comparable  with  w  or  the  spectrum  spread,  then  the  Fokker-Pl  anck 
diffusion  approximation  breaks  down.  This  begins  to  happen  in  N2  at 
average  energies  of,  say,  1  or  2  eV  when  energetic  tail  electrons  can 
excite  electronic  transitions  and  w^  -  7  eV,  or  at  -  .01  eV  where  rota¬ 
tional  jumps  have  comparable  energy.  There  is  an  interesting  intermediate 
regime  where  the  dominant  w^  is  sufficiently  smaller  than  w  but  not  negli¬ 
gible.  This  appears  to  be  the  case  in  n2  for  E/N 0  >  .01  Td,  giving  the 
swarm  energies  greater  than  rotational  jumps,  and  where  vibrational  tran¬ 
sitions  dominate.  Even  at  1  or  2  eV,  electronic  transitions  induced  by 
the  tail  are  only  a  part  of  the  energy  loss,  and  the  diffusion  approxima¬ 
tion  is  still  not  bad. 

One  important  feature  of  keeping  M  is  that  this  term  allows  some 
electrons  to  lose  energy  slower  than  the  mean,  thereby  supporting  the  high 
energy  tail  of  the  spectrum,  as  seen  graphically  in  Figure  lb.  This  will 
greatly  increase  the  avalanche  rate  over  the  CSDA. 

We  have  still  neglected  third  and  higher  derivatives  in  the 
expansion  Eq.  (22).  The  main  properties  of  any  normalized  distribution 
with  a  decaying  tail,  such  as  F  or  g  ■  /w  F,  are  its  average  energy  and 
the  “width"  of  the  function  about  the  average.  The  first  derivative  term 
L  determines  the  average  energy,  and  the  second  derivative  term  M  deter¬ 
mines  the  width,  together  with  corresponding  terms  from  the  field  contri¬ 
bution.  Higher  neglected  terms  could  affect  only  the  detailed  shape  of 
the  function.  For  example,  the  third  derivative  would  alter  the  asymmetry 
of  the  spectrin  about  the  mean,  and  the  fourth  derivative  would  set  the 
kurtosis.  The  transport  coefficients  of  interest  are  weighted  integrals 
over  the  spectrum,  and  the  details  of  the  spectrin  tend  to  integrate  out. 
Thus  there  is  physical  reason  to  expect  that  stopping  with  L  is  insuffi¬ 
cient  (CSDA)  but  stopping  with  M  is  sufficient  to  calculate  a  spectrum 


which  reproduces  bulk  transport  coefficients  with  some  acceptable  accur¬ 
acy.  These  arguments  may  be  less  true  for  the  avalanche  coefficient  at 
low  temperatures  which  requires  an  accurate  high  energy  tail. 

CORRECTIONS  FOR  GAS  HEATING 

Equations  (23)  for  the  inelastic  parts  of  L  and  M  are  expres¬ 
sions  for  the  gas  at  rest.  Thermal  effects  were  mentioned  only  in  connec¬ 
tion  with  the  elastic  part  of  L  in  Eq.  (12).  In  fact  the  Davydov  correc¬ 
tion  for  thermal  motion  of  molecules  shows  that  collision  calculations  are 
to  be  modified  by  replacing  g  by: 

g  ♦  q+T/wL  f_El  (30) 

9  3w  /w 

or,  in  terms  of  F(v): 

F  *  F  +  T  l£  (31) 

q  3  w 

Therefore  thermal  qas  motion  is  taken  into  account  by  replacing  the  flux  $ 
in  Eq.  (28)  by 

$  =  &2  E2  fg  -  2.9_  (wg)  1  -  N0vL[g  +  T  /w  -2—  f— )  1 

mv  3  3w  9  3w  /w 

m 

(32) 

-  JL_  fNnvMig  +  T  /w  2_  f—  111 
3w  9  aw  /w 

The  steady  state  spectrum  is  that  for  which  $  =  0.  When  E  +  0  as  well, 
the  solution  of  (32)  clearly  reduces  to  a  Maxwellian,  as  it  must,  for  then 
g  +  T  /w  d(g//w)/dw  =  0. 


SECTION  4 

EQUATION  FOR  TIME  DEPENDENCE  OF  ENERGY  SPECTRUM 


* 


» 


* 


Dropping  the  Davydov  terms,  equations  (27)  and  (28)  (or  (32) 
with  =  0)  control  the  time  evolution  of  the  energy  spectrum  in  the  ab¬ 
sence  of  sources  and  sinks: 

(e_2E 2  Fg  -  2  \(wg]T  -  N.vLg  -  N  .  D(vgM)t  =  Q  (33) 

>t  ^w  rnvm  '  3  *  0  3  *  1 

m 

Sources  and  sinks  are  to  be  added  to  the  right  hand  side  for  a  specific 

-V 

problem  as  appropriate.  Eq.  (33)  is  valid  for  E  varying  slowly  compared 
with  the  momentum  transfer  freauency  but  on  any  time  scale  relative  to 
the  mean  energy  transfer  frequency.  Since  j  >  lO'Vsec  in  sea  level  air, 
Eq.  (33)  will  apply  to  many  problems  of  interest.  For  example,  following 
a  pulsed  ionization  source,  for  which  the  starting  spectrum  is  the  initial 
condition  to  Eq.  (33),  this  equation  determines  the  relaxation  of  the 
spectrum  toward  equilibrium  in  the  electric  field  and  warm  gas. 

For  situations  where  E  =  0,  the  resultinq  equation 

-  N0  I vL g  +  ■  =  0  (34) 

i  c  }w  v 


determines  the  relaxation  of  q  from  an  initial  condition  toward  thermal 
equilibrium.  Without  the  E  term,  Eq.  (34)  is  no  lonqer  limited  by  v^.^ 

<<  v . ,  ,  (two  term  Leqendre  series),  and  so  applies  to  electrons  of  any 

thermal 

energy.  It  includes  mean  enerqv  loss,  straqqlinq,  and  thermal  heatinq 
(through  thermally  populated  rotational  states  in  l  and  M) . 
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If  straggling  were  to  be  neglected  (CSDA),  the  resulting  equa¬ 
tion 

#  -  N0  ( vgL )  =  0  (35) 

5t  Tw 

would  show  simple  convection  in  energy  space,  and  laminar  flow.  That  is, 
electrons  of  higher  energy  would  never  overtake  ones  of  lower  energy. 
Straggling  destroys  laminar  flow,  and  permits  passing  in  energy  space,  as 
actually  occurs  in  a  random  walk,  and  spreads  out  the  spectrum. 


SECTION  5 

STEADY  STATE  SPECTRUM 


In  steady  state,  the  solution  for  F  -  (m/v)g  of  Eq.  (33)  with 
3;g/  >t  =  0  is 


F 


F  :■  exp  j 


w  L  +  I  (wM) 
w  Tw 

o  ( eE/N  q) 2  +  y 


dw  • 


(36) 


where  Fo  is  a  normal i zat ion  constant  set  by  Fq.  (4).  This  spectrum  is  an 
improvement  over  the  CSDA 


3  7  L 

F  =  Fj  expi  -  f  m  dw’>  (37) 

(eE/No)2 

which  results  when  M  is  neglected.  Both  have  a  auite  simple  functional 
form. 


THE  DRUYVESTEYN  SPECTRUM 


The  Druyvesteyn  spectrum  is  recovered  from  Eq.  (37)  in  the 

special  case  of  j  =  constant  and  elastic  collisions  only, 
m 


L 


L  =  2m  w 

el  ^  m 


(38) 


m 
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Then,  F  reduces  to 


* 


.  3m  of  w  , 

F  =  F  o  ex  pi-  - 2 - 4  (39) 

M0(eE/N0)2 

In  this  same  special  case  the  more  correct  Eq.  (36)  is  essen¬ 
tially  the  same  as  (39)  because,  for  elastic  collisions  only,  M  -  (2m/ 
M0)wL  <<  wL,  and  the  M  correction  is  negligible.  Thus  the  improved  spec¬ 
trin  does  not  mar  already  good  agreement  with  data. 

It  is  when  inelastic  collisions  are  important  that  the  differ¬ 
ence  between  (36)  and  (37)  is  significant. 

THE  NEW  TERMS 


In  Section  6  it  will  be  shown  that  the  two  terms  (eE/N0)2/3  am 
and  M  in  the  denominator  of  Eq.  (36)  can  be  identified  as  diffusion  coef¬ 
ficients  in  energy  space  due  to  elastic  scattering  in  the  electric  field 
(first  term)  and  due  to  spread  in  energy  loss  about  the  mean  (M  term). 

The  total  diffusion  coefficient  is  just  the  sin  of  the  two. 

The  first  term  L  in  the  numerator  is  the  mean  rate  of  energy 

loss  (or  gain,  if  L  is  negative)  due  to  collisions.  The  second  term 

(l/w)d(wM)/dw  is  a  mean  loss  (or  gain)  rate  due  to  energy  variation  of  the 
diffusion  coefficient  M. 

At  energies  below  vibrational  excitation,  and  neglecting  elastic 
scattering,  only  rotational  states  contribute  to  L  and  M.  We  shall  use 
the  Gerj uoy-Stein  rotational  cross  sections  to  evaluate  this  contribution 
(Appendix  A).  States  Aj  =  2  apart  are  spaced  by  ( 4J-2 ) B 0 ,  where  B0  is 

the  rotational  constant.  The  cross  section  for  de-excitation,  a 
diverges  as  /4JB0/w  when  w  <  4JB0.  In  nitrogen,  B0  =  2.5*10-4 


J,  J-2 
eV.  This 
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causes  a  formal  divergence  in  the  numerator  term  (l/w)d( wM)/dw,  and  effec¬ 
tively  makes  the  soectrum  unreliable  for  w  <  10"  3  eV,  and  the  calculated 
transport  coefficients  inaccurate  for  E/N^  <  10“  2  T  .  When  w  <  4JB  :  dif¬ 
fusion  is  no  longer  a  good  approximation.  For  numerical  calculations  we 
have  simply  smoothed  out  the  divergence  by  multiplying  this  numerator  term 
by  w/(w+Bq).  The  inaccuracies  in  characteristic  energy  and  drift  velocity 
will  be  seen  in  Section  7. 

THE  LOW  ENERGY  THERMAL  LIMIT 


Eq. 


For  vanishingly  small  electric  fields,  the  remaining  terms  in 

(36)  should  make  F  reduce  to  a  Maxwellian  at  T  .  Thus  we  should  have 

9 


L  +  w-1  T(WM)/*w  =  !_ 
M  T9 


(40) 


or,  solving  for  M, 


=  -(w'-w)/Tn 

M  =  1  j  e  w'l(w' )dw'  (41) 

w  w 


M  must  be  a  weighted  integral  of  L  over  an  energy  span  about  T  . 

Since  T  enters  L  and  M  only  through  rotational  populations  (we 

g 

neglect  vibrational  temperature),  Eq.  (40)  should  hold  if  we  consider  only 
the  rotational  contributions.  Upon  constructing  L  and  M  (Appendix  A)  one 
finds  Eq.  (40)  is  in  general  not  true.  Rotational  levels  with  J  =  10  or 
larger  exist,  and  so  energy  jumps  4JB  i  =  4 *10 *2 . 5  *10" 4  =  .01  eV  in  N: 
occur,  encroaching  on  electron  thermal  energies  ~  .02  eV.  That  is,  the 
diffusion  approximation  breaks  down  at  these  low  energies,  and  Eq.  (36)  is 
not  a  good  approximation  to  the  spectrum. 
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However  in  a  hypothetical  gas  with  Bj  +0,  all  energy  jumps  are 

small,  and  Fa.  (36)  is  correct.  For  such  a  qas,  Eq.  (40)  should  hold.  In 

Appendix  A,  using  the  Gerjuoy-Stein  rotational  cross  sections,  we  take  the 
limiting  expressions  of  L  and  M  for  B0  *  0,  and  show  explicitly  that  Equa¬ 
tions  (40)  and  (41)  are  indeed  true,  independent  of  all  cross  section 
parameters.  According  to  equilibrium  statistical  mechanics,  this  "must" 
happen,  but  it  is  comforting  to  see  it  follow  from  a  direct  calculation 
using  actual  cross  sections. 

Using  the  complete  L  and  M  for  N2,  the  actual  spectrum  of  Eq. 

(36)  with  E  =  0  is  shown  in  Figure  2.  Above  about  .01  eV  it  is  indeed 

Maxwellian  at  the  correct  temperature,  300  K.  (The  departure  in  the  high 
energy  tail  above  0.3  eV  is  because  we  assume  the  gas  is  in  the  vibration¬ 
al  ground  state).  But  near  .01  eV  there  is  a  non-Maxwel  1  i an  dip  due  to 
"large"  energy  jumps  not  well  modelled  by  the  Fokker-Pl anck  form.  This 
irregular  behavior  goes  away  as  B0  -*•  0. 


ure  2. 


Spectrum  in  N2  from  Equation  (36)  when  E=0.  Dip  near  .01  eV 
is  due  to  rotational  transitions  of  comparable  energy.  Gas 
at  300K . 


SECTION  6 

RELATION  TO  ORDINARY  DIFFUSION  THEORY 


The  functional  form,  Eq.  (36),  for  F  has  a  simple  interpretation 
in  terms  of  elementary  diffusion  and  convection  theory. 

Let  a  quantity  ip  diffuse  and  convect  in  one  dimension  x  with 
diffusion  coefficient  D  and  convection  velocity  V: 


11  +  i_  ( v t|>)  =  JL  D  14  (42) 

at  ax  ax  ax 

We  shall  presently  identify  x  with  w  and  ^  with  g.  V  and  D  may  be 
functions  of  x.  In  steady  state,  the  elementary  diffusion-convection 
equation  (42)  shows  that  the  flux  must  be  constant  in  t  and  x: 

VP  -  D  11  =  constant  =  0  .  (43) 

ax 

We  assume  '>  and  its  derivative  vanish  at  some  end  point,  so  the  constant 
is  taken  to  be  zero.  Then 

1  11  =  1  .  (44) 

ip  3x  D 


In  steady  state  the  (prescribed)  convection  velocity  and  diffusion 
coefficient  determine  the  logarithmic  derivative  of  the  dependent  variable 
by  Equation  (44).  Consequently  is  distributed  according  to 


where  *  is  a  normal  izat' 


V 


I 

V  2 


-  1-  (N 0vM) 
3w 


=  energy  drift  due  to  aui 

3w 


(52) 


Diffusion  coefficients  have  dimensions  energy2/t ime,  and  velocities  have 
dimensions  of  energy/time.  The  energy  dependence  of  and  D-  cause  a 
mean  flux  (V2  and  v])  because  of  the  way  they  appear  in  Eq.  (46)  under 
the  32/3w2  operator.  Eq .  (46)  can  then  be  written 


wi  th 


and 


la  +  J_  (Vg)  =  i_  (D  ia) 

3t  3w  3w  3w 


V=vj+vl+v{+v] 

=  N0 v  rl£g/.Nol2  -  L  -  —  —  (vM)  1 

3mv2a  v  5w 

m 


D  =  De  +  Di 


=  n0v  +  Ml 

3  a 


(53) 


(54) 


(55) 


Following  Eq.  (45)  the  steady  state  solution  of  Eq.  (53)  is 

g  =  g0  exp  ( f  1  dw)  (56) 

Forming  the  ratio  V/D  from  Equations  (54)  and  (55),  we  add  and  subtract 
M/2w  in  the  numerator,  and  use  v" 1 ( 3/3w) ( vM)  +  M/2w  =  w" 1  (  3/ 3w)  ( wM)  to 
f  i  nd 


/ 
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(57) 


* 


L  +  I  i_  (wM) 

V  =  1_  .  w  ?w 
D  2w  E^f.  +  m 

3  a 

m 

where  E'  =  eE/N0  is  the  reduced  electric  field.  Upon  integrating  and 
exponentiating,  the  first  term  (l/2w)  becomes  /w  ,  and  the  remaining 
factor  is  therefore  F, 

F  =  Fn  expl-r  L  W1  9(wM)/3w  riwj  (58) 

0  E  '  2  /  3  a  +  M 
u  m 

as  in  Eq.  (36).  Thus  the  functional  form  for  the  steady  state  energy 
spectrum  is  an  elementary  consequence  of  ordinary  convection  and  diffusion 
in  energy  space,  as  is  the  time  dependent  equation  (46)  or  (33). 

This  can  be  seen  as  well  by  writing  the  flux,  Eq.  (28),  in  terms 
of  F.  After  some  algebra  it  becomes 

*  =  -  ^0  f rL  +  1  1  F  +  f—  +  Ml  £.1  ,  (59) 

m2  w  3w  3am  3w 

the  coefficients  of  F  and  8F/9w  being  the  convection  velocity  numerator 
and  diffusion  coefficient  denominator,  respectively,  of  the  integral  in 
Eq.  (58). 
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SECTION  7 

THE  CASE  OF  NITROGEN 


We  apply  these  ideas  to  N2.  L  and  M  are  constructed  from  cross 
section  data.  The  steady  state  spectrum,  Eq.  (36),  and  swarm  parameters 
are  then  computed. 

CROSS  SECTIONS,  AND  L  AND  M  FUNCTIONS 

Electron-molecule  cross  sections  are  compiled  in  a  number  of 
places  (refs  9  -  12).  Reference  13  is  a  recent  critique  of  rotational  and 
vibrational  excitations.  Here  we  use  the  Ger juoy-Stein  expressions  (ref 
14)  for  rotational  cross  sections,  with  the  N2  quadrupo le  moment  Q  =  1.04 

in  units  of  ea2. 

0 

For  the  vibrational  and  electronic  cross  sections,  we  use  a  set 
compiled  by  Archer  (ref  15)  that  has  been  in  use  in  the  DNA  Auroral 
physics  community  for  years.  These  vibrational  excitations  are  based  on 
the  measurements  by  Schulz  as  normalized  by  Haas.  A  low  energy  (0.29  - 
1.7  eV)  tail  for  the  v  =  0  -  1  transition  was  added  as  suqgested  by 
Englehardt,  Phelps,  and  Risk  (Ref.  6).  Numerical  values  for  the  tail  were 
taken  from  the  second  table  on  page  64  of  Kieffer  (Ref.  16).  The  elec¬ 
tronic  cross  sections  agree  quite  well  with  those  tabulated  by  Wadzinski 

and  Jasperse  ( ref  11 ) . 

L  and  M  for  10" 3  <  w  <  10 2  eV  are  shown  in  Figure  3a.  L  becomes 

negative  for  w  <  Tq  ~  .025  eV,  because  a  lower  energy  electron  gains 


energy  on  the  average  by  collisions  with  the  gas  at  this  temperature  from 
super-el ast ic  rotational  transitions.  M,  of  course,  is  always  positive. 
Over  10- 3  to  10  eV,  these  functions  vary  by  more  than  six  orders  of  magni¬ 
tude.  The  ratio  M/L  is  shown  in  Figure  3b.  It  is  much  more  slowly  vary¬ 
ing  than  M  or  L  itself. 

THE  ENERGY  SPECTRUM 


Figures  4  and  5  show  the  energy  spectrum  F,  Eq(36),  for  E/No  = 
10  Td  and  100  Td,  respectively.  1  Td  =  10" L/  V-cm2.  The  high  energy  tail 
is  markedly  higher  than  in  the  CSOA  which  is  also  shown.  Any  process  with 
significant  cross  section  above  the  CSDA  tail  will  have  markedly  different 
rates  with  the  two  spectra. 

Figure  5  also  shows  a  spectrum  from  the  calculation  by  Pitchford 
and  Phelps  (Ref.  17),  involving  a  numerical  solution  of  the  Boltzmann 
Equation  with  complete  collision  integral.  These  authors  used  a  different 
cross  section  set,  so  the  comparison  is  not  meaningful  in  a  precise 
sense.  Nevertheless  the  differences  are  noteworthy.  In  particular,  we 
overestimate  the  high  energy  tail  and  underestimate  the  low  energy 
population  below  2  eV.  This  is  due  to  the  local  diffusion  approximation. 
Electronic  transitions  allow  high  energy  electrons  to  jump  to  very  low 
energies  more  quickly  than  in  the  diffusion  model.  This  excessive  tail 
will  also  cause  a  large  avalanche  coefficient  discussed  later. 

Eq.  (36)  also  shows  a  curious  broad  plateau  around  3.5  to  7  eV, 
absent  in  the  more  complete  calculation.  This  too  is  likely  due  to 
electronic  transitions  which  permit  electrons  to  jump  over  this  energy 
range.  In  general  it  is  clear  the  average  swarm  energy  at  100  Td  will  be 
computed  to  be  too  large. 
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TRANSPORT  COEFFICIENTS 


The  open  circles  in  Figure  6  show  the  drift  velocity  calculated 
with  E q ( 36 )  superimposed  on  data  compiled  by  Dutton  (ref  8).  It  is 
slightly  too  small  below,  say,  3*10“ V- cm2.  Here  the  Fokker-P 1 anck 
approximation  is  in  jeopardy  because  rotational  transitions  become 
comparable  to  swarm  energies.  Refer  to  the  spectrum  in  Figure  2  at  zero 
field.  Average  energies  are  slightly  too  large  as  seen  momentarily. 

The  drift  velocity  is  also  too  small  near  10" lo  or  10"  V- 
cm4.  Here  the  spectrum  is  similar  to  that  in  Figure  5,  and  averaqe  ener¬ 
gies  are  again  too  high. 

At  stronger  fields  the  drift  velocity  returns  to  good  agree¬ 
ment.  Rut  at  these  larger  energies  the  mono-energetic  momentum  transfer 
collision  frequency  in  N2,  N0v cr^  ,  is  slowly  varying  with  energy,  and  the 
drift  velocity  is  only  weakly  dependent  on  the  spectrum. 

At  fields  less  than  0.1  Td,  we  find  drift  velocities  computed  in 
the  CSOA  to  be  more  than  20%  too  small,  and  characteristic  energies  more 
than  30%  too  large.  But  in  stronger  fields,  we  find  CSDA  drift  velocities 
to  be  more  like  5  to  12%  too  small,  and  characteristic  energies  5  to  13% 
1 arger  than  in  the  diffusion  approximation. 

The  swarm  characteristic  energy  is  shown  in  Figure  7.  It  is 
slightly  too  large  below  about  3*10-i:l  V-cm",  and  is  significantly  too 
large  above  10_  Io  V-cm2,  accounting  for  the  lowered  drift  velocities.  As 
mentioned,  this  is  due  to  finite  rotational  and  electronic  energy  jumps  in 
these  two  regions. 


The  swarm  momentum  transfer  ( v  )  and  energy  transfer  ( vu )  colli¬ 
sion  frequencies  are  shown  in  Figure  8,  superimposed  on  a  graph  from  Frost 


Elastic  (krn/N)  and  Energy  Exchonge  (14/N) 
Collision  Frequency  (cm3/sec) 


Figure  8 


Swarm  momentum  transfer  (v  /N)  and  energy  transfer 
( v  /N)  collision  frequencies.  The  large  open  circles 
are  our  calculated  values,  superimposed  on  a  graph  from  Frost 
and  Phelps  (Ref.  18).  The  solid  curves  are  averaged  experi¬ 
mental  data  at  300  K.  The  other  four  sets  of  points  are 
various  theories  within  Ref.  18. 


and  Phelps  (Ref.  18).  They  are  defined  by 


V  = 


eE 


mv 


drift 


eE  v 


vu 


drift 


D/  u  -  T 


(60) 


(61) 


where  u  is  the  mobility,  and  D/p  is  the  swarm  characteristic  energy. 
(Ref.  18  uses  for  the  swarm  frequency,  a  symbol  we  have  reserved 

for  the  mono-energetic  collision  frequency.)  The  agreement  is  striking 
except  for  near  1  eV.  The  dependence  of  on  and  0/ p  accen¬ 

tuates  the  dependence  on  the  swarm  average  energy. 


The  non-self-consistent  avalanche  coefficient,  a(cm-1)>  was  com¬ 


puted  by 


m  = 


1 


'drift 


/dwgv^. 


(62) 


where  ^  is  the  total 


N2  ionization  cross  section.  It  is  shown  as 
open  circles  in  Figure  9,  superimposed  on  data  of  Dutton.  Below  about  300 
Td  the  CSDA,  shown  as  crosses,  is,  of  course,  woefully  low.  Eo.  (36)  is 
too  high  by  more  than  an  order  of  magnitude  at  100  Td,  a  discomforting 
amount.  The  spectral  tail  seen  in  Figure  5  is  responsible. 

DISCUSSION 

With  some  exceptions  (avalanche  rate,  characteristic  energy  when 
near  several  eV)  the  Fokker-Pl anck  approximation  to  the  energy  spectrum 
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appears  to  be  able  to  reproduce  transport  coefficients  in  N2  with  reason¬ 
able  accuracy.  Other  gases  are  under  study.  Inaccuracies  show  up  when 
molecular  energy  jumps  are  not  small  compared  with  swarm  energies,  pre¬ 
cisely  when  the  diffusion  approximation  becomes  unreliable.  These  fea¬ 
tures  are  understandable  when  it  is  realized  that  it  is  the  spread  of 
energy  gain  and  loss  about  the  mean  rates  that  actually  determine  the 
spectrum. 

At  present  the  ultimate  utility  of  the  diffusion  approx imat ion 
is  unclear.  It  depends  in  part  on  the  success  in  other  qases.  It  may  be 
quite  useful  for  rough  calculations,  it  being  possible  to  obtain  spectra 
and  all  transport  coefficients  merely  by  numerical  quadratures.  The  time 
dependent  equation  (33)  may  be  especially  useful,  being  a  significant 
improvement  over  the  CSOA,  but  not  substantially  more  complicated. 


APPENDIX  A 

ROTATIONAL  CONTRIBUTION  TO  L  AND  M 


The  rotational  contribution  to  ,  Eq.  (21),  is 


Crot  '  /=(,  -  (9V»J,J«>w1 


(A-l) 


*  "  Njr<’,»J,J-2)w-»,  ,  ,  -  (qVc,J,J-2>w1 


J,  J-2 


J,  J-2 ' w 


since  only  J±2  are  accessible  from  J  for  a  quadrupole  interaction.  The 
second  line  in  Eq.  (A-l)  is  collisions  of  the  second  kind,  in  which  an 
electron  gains  enerqy.  Nj  is  the  number  of  molecules  in  the  J-th  state. 


NJ  =  No 


pj 


pj  * 


Pj  =  3(l+a) (2J+1 )  , 


Ej  =  J(J+1 )B0  , 


(A-2) 


( A-3) 


and  a=0  for  odd  J,  a=l  for  even  J.  The  transition  energies  are 


JJ,J+2  =  E.J+2  ‘  EJ  =  (4J+6>Bc 


( A-4 ) 


WJ , J-2  "  EJ  "  E J-2  =  ^4J"2^Bf 


The  cross  sections  for  excitation  and  de-excitation  are  taken  as  those  of 
Gerjuoy  and  Stein  (Ref.  14). 


where 


,,  =  u±au±n  Oou  -  (4^)Bqi1/2 

J»J+2  (2J+3) (20+1)  0  w 


O,  ,  -  •  11ML  „,ri»  <4J-g)B0]1/? 

J’J-2  (2J-1 ) (2J+1 )  0  w 


a0  =  8-nQ2  a^/15 


( A-5 ) 

(A-6) 


( A- 7 ) 


a0  =  Bohr  radius  =  .529xl0-8  cm 


and  Q  is  the  molecule  quadrupole  moment  in  units  of  ea2.  Q=1.04  for  N; 
B0  is  the  rotational  constant,  =  2.5xl0"4  eV  in  N2. 


The  expressions  for  L  and  M,  Equations  (23)  are 


-rot  =  V  NJwJ,J+2°J,J+2  "  V  NJWJ, J-2°J, J-2 


(A-8) 


M  =  i  V  N.w2  .  „a  •>  ,  „  +  i  V  N.w2,  ,  0cj,  ,  - 
rot  2  J  J,J+2  J,J+2  2  J  J,J-2  J,J-2 

0  2 


(A-9) 


At  300’  it  takes  20  or  30  terms  for  the  sums  to  converqe  (by  which  point 
E j/T g  =  J(J+1 )B0/Tg  >  5). 


Below  the  threshold  for  exciting  the  lowest  state,  J=0  to  J=2,  w 
<  6B0,  L  is  obviously  negative.  The  electron  can  only  qain  energy.  All 
terms  together  actually  make  L  .  negative  for  w  <  T  . 


THE  SMALL  B0  LIMIT 


We  now  obtain  the  limiting  expressions  for  l_rot  and  M  t  when 


B0  <<  T  ,  and  B0  <<  w.  Let 


a  =  —  <<  1  ,  y  =  «  1  . 

Tg 


( A-10) 


Inserting  the  variables  in  Eq.  (A-8)  and  cancelling  common  factors  in  the 
numerator  and  denominator  we  have 


L  =  2Acf0 B0  y  (l+a)e'J(J+1)a  (J+2) (J+l )/l-(4J+6)y 


oo 

-  2Aa0B0  y  (l+a)e'j(j+1)ct  J(J-1)/1+  (4J-2)’ 


( A- 11) 


where  A  =  3N0 /i jP jexp(-J( J+l)o)  is  a  normalization  constant. 

Many  terms,  up  to  J  -  /5/a,  will  contribute.  Take  the  average 
value  of  (1+a)  to  be  3/2,  and  keep  first  order  in  y : 


L  =  3A0OBO  y  (J+2)(J+l)e~J(J+1)ct  f  1- (2 J+3 ) y ] 


-  3Ao0B0  y  J(J-l)  e"J(J+1)a  ri+(20-l)yl 


( A- 12) 


The  sums  will  behave  as  inverse  powers  of  a,  so  in  the  first  sum  the  J=n 

oo 

and  J=1  terms  are  negliaible  compared  with  the  remaining  i  .  Combining 

.  .  2 
the  two  remaining  sums  I  , 


_L -  =  v  i;ai  (J+2)(J+1)  -  J(J-l) 

3Aa0B0  2 

-f (J+2)(J+l)(2J+3)  +  J(J-1)(2J-1)  V 

=  2  V  e-J (J+1 )a{ 2j+i  .  r 2J 3+3J 2+7J+3 ly 1 
2 


Higher  powers  of  J  lead  to  higher  powers  of  cT1//2  in  the  sums, 
leading  terms, 


1 _  =  4  ”  e'J(J+1)arJ-^ 


1 


3Aan  B 


0U0 


The  sums  can  be  approximated  by  integrals 

V  Je‘J(J+1'a  -  r  xe'aX  dx  -  f  xe'aX  dx  = 


1 

2a 


"  j3e-J(J+l)a 


3o-^>X  Hv  =  1 


f  xse  dx  = 


2a‘ 


Then 


_ k _ =  2.  n  _  x)  =  a  -  lai 

3ActqBq  a  a  B0  w 


The  sum  in  A  is 


V  p  e-J(J+l)a  „ 

J  J  2a 


SO 


A  =  3N0/EJPJe"J(J+1)a  =  |  aN o 

Equations  (A-16)  and  (A-17)  then  show  the  limiting  form  for  L  as 
be 


T 

L  =  4N0a0B0fl  -  -£) 


( A- 13) 

Keeping 

( A- 14 ) 

(A-15a) 

(A-lAb) 

(A-16) 

(A-17) 

Bg  0  to 


( A-18) 


Electrons  with  w  <  Tq  qain  energy. 


A  parallel  calculation  for  M  shows 


M  =  4NnanB«T 


( A- 19) 


Equations  (A-18)  and  ( A- 19)  are  explicit  expressions  for  L  and  M  as  B0  * 
0.  In  this  same  limit  the  diffusion  approximation  to  the  spectrum  should 
be  val id . 

REDUCTION  OF  SPECTRUM  TO  MAXWELLIAN 

The  steady  state  spectrum  ‘Eq.  (36)  must  become  Maxwellian  when 
E  =  0.  This  led  to  the  requirement  (40)  or  (41).  We  now  show  these 
requirements  are  indeed  met  when  B0  +  0  for  then  energy  jumps  are  small 
and  the  diffusion  expressions  are  good. 

With  M  from  Eq.  (A-19),  and  L  from  (A-18),  find 


1  (wM)  =  M 
w  ?w  w 


(A-20) 


L  =  M_  -  M 

T  w 

q 


( a-21 ; 


L  +  w~l  ?(wM)/3w  =  1_ 
M  T 


( A-22 ; 


as  required.  Likewise  the  integral  relationship  (41)  between  L  and  M  is 
easi ly  shown  to  hold. 


This  shows  by  explicit  calculation  that  for  small  enerqy  jumps 
the  diffusion  approximation  reduces  correctly. 
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